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The chiral-odd generalized parton distribution (GPD), or transversity GPD, of the nu- 
cleon can be accessed experimentally through the photo- or electroproduction of two 
vector mesons on a polarized nucleon target, j^N — > pip 2 N', where p\ is produced at 
large transverse momentum, p 2 is transversely polarized, and the mesons are separated 
by a large rapidity gap. We predict the cross section for this process for both transverse 
and longitudinal p 2 production. To this end we propose a model for the transversity GPD 
HT{x,l;,t), and give an estimate of the relative sizes of the transverse and longitudinal 
p 2 cross sections. We show that a dedicated experiment at high energy should be able to 
measure the transversity content of the proton. 
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1 Introduction 



Accessing the chiral-odd quark content of the proton has a long history [1] which has 
been detailed elsewhere [2]. The particular case of the chiral-odd generalized parton 
distributions is interesting in itself [3] . Following the previous work in Ref . [4] , we consider 
here the process 3 

(*) o + /i \ 

Tl/tP PlPl/t h I 1 ) 

shown in Fig. ^ that is, virtual or real photoproduction on a proton p, which leads 
via two-gluon exchange to the production of a longitudinally polarized vector meson p° L 
separated by a large rapidity gap from another longitudinally or transversely polarized 
vector meson p^, T and the scattered neutron n. We consider the kinematical region 
where the rapidity gap between p + and n is much smaller than the one between p° and 
p + , that is the energy of the system (p + — n) is smaller than the energy of the system 
(p° — p + ) but, to justify our approach, still larger than baryonic resonance masses. 

We have previously shown [4] that in such kinematical circumstances the Born term 
for this process is calculable consistently using the collinear QCD factorization method. 
The final result is represented as an integral (over the longitudinal momentum fractions 
of the quarks) of the product of two amplitudes. The first one describes the transition 
7M — > p° L via two-gluon exchange which can be also viewed as the Born approximation of 
a hard Pomeron. The second one describes the two gluon (Pomeron)-proton subprocess 

— > p + n which is closely related to the electroproduction process 7* iV — > pN' where 
collinear factorization theorems allow separating the long distance dynamics expressed 
through the GPDs from a perturbatively calculable coefficient function. The hard scale 
appearing in the process in Fig. ^ is supplied by the relatively large momentum transfer 
in the two-gluon channel, i.e. by the virtuality of the "Pomeron" . 

The first process we will calculate is the one with both vector mesons longitudinally 
polarized, 

7 W (?) P(P2) - P°M pUp p ) n(p' 2 ) , (2) 

which involves the emission of two gluons in the 7 — > pl transition . We choose a charged 
vector meson p + to select quark antiquark exchange with the nucleon line. The second 
process we are interested in is the one involving the chiral-odd GPD, e.g. 

T W(g) p(j)2 ) pl ( ?p ) p+( Pp ) n ( p ' 2 ) , (3) 
which is the main motivation for the study of this two meson production process. 

2 Kinematics 

Let us first summarize the details of the kinematics of the process, restricting to real 
photoproduction. We introduce two light-like Sudakov vectors p\ and P = 1/2 (p 2 +P2>)- 

3 The process j*p — > p\p\y may equally well be discussed along the same lines, since its amplitude is 
described by exactly the same graphs, thanks to charge conjugation invariance which forbids contribution 
of the gluonic GPDs. 
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Figure 1: Factorization of the process ^Jt P ^ P°l pX/t n in the as y mme tric kinematics 
discussed in the text. P describes the two gluon exchange (Born approximation of the hard 
pomeron). 



We also introduce the auxiliary variable S = 2pi P, related to the total center-of-mass 
energy squared of the 7*p— system, s = (q + P2) 2 , as s + Q 2 = (1 + £)S. The momenta 
are parametrized as follows : 

O 2 

1 ao 

P^afi + ^P^-pfr, a = l-a, 
ao 

^ = (1+0^, ^, = (1-0^, (4) 

where Q 2 = —q 2 is the photon virtuality, and £ is the skewedness parameter which can 
be written in terms of the invariant mass si of the two mesons as 

Si + Q 2 2 p 2 

£ = 9g » s 1 = (q p +p p ) = —. (5) 

10 aa 

The p + -meson-target invariant mass equals 

S2 = {p p + p 2 >Y = Sa (1-0 . (6) 

The kinematical limit with a large rapidity gap between the two mesons in the final state 
is obtained by demanding that s\ be very large, of the order of S ~ s, 

Sl = 2S£, Sl »p 2 , (7) 

whereas S2 is kept of the order of p 2 , 

f(l-0, (8) 
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large enough to justify the use of perturbation theory in the collinear subprocess ¥p — > 
p + n and the application of the GPD framework. 

In terms of the longitudinal fraction a the limit with a large rapidity gap corresponds 
to taking the limits 

a — > 1 , asi — > p 2 , £ ~ 1 . (9) 

We choose the kinematics so that the nucleon gets no transverse momentum in the 
process. The t-dependence of the GPDs is not known in detail, but the sum rules which 
relate it to the nucleon form factors imply a strong decrease of the cross section with 
increasing —t. Thus, taking t = i m j n yields simpler formulas and amounts to most of the 
cross section. One may however allow a finite momentum transfer, small with respect 
to \p\, correcting expressions of the amplitude with contributions containing additional 
GPDs. 

Let us repeat that the role of the main hard scale in the processes under discussion 
below is played by the virtuality p 2 = —p 2 , which is the large momentum transfer in the 
two-gluon exchange channel. 

The 3-particle phase space needed to obtain the differential cross section is computed 
through the standard recurrence relation, which gives 

dPS 3 {s,p 2 ,p p ,p 2 >) = 7^7— a dp 2 T dtds 2 d§ , (10) 

where $ is the angle between the (jp°) plane and the (p, n, p + ) plane. At t — t min , this 
angle is irrelevant and the phase space factor simplifies to 

dPS 3 (s,p 2 ,Pp,p2>) = 10Q 3 \ n -rdp 2 T dtd£ . (11) 

lzo7r° s 4 (1 — 4 > 

The differential cross section at the minimal value of —t = —t m i n then reads 

do 1 
dp\dtdi = 256vr 3 £(l -4> 2 '" M ' ' ^ 

3 The scattering amplitude 

The scattering amplitude M. of the process (j2J) using the standard collinear QCD factor- 
ization method is written in a form suggested by Fig. Q a s ; 

i i i 

Ai ~ J dz J du J dxTfl(xi,u, z) H q (x,£,0)(f)p+(u)(f)po(z) . (13) 

q=u,d oo-l 

Here H g (x,l;,0) is the generalized parton distribution of parton q in the target at zero 
momentum transfer and x + £ and £ — x are (see Fig. [IJ the momentum fractions of the 
quark and antiquark emitted by the target (since it turns out that our kinematics selects 
x < £, the second parton is interpreted as an emitted antiquark). 4> p +(u) and (fi p o(z) are 
the distribution amplitudes (DA) of the p + -meson and p°-meson, respectively, and u and 
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z are the corresponding lightcone momentum fractions of the quark in the meson. We 
will also use the shorthand notation z = 1 — z and u = 1 — u in the following. T^(x, u, z) 
is the hard scattering amplitude (the coefficient function). For clarity of notation we 
omit in Eq. (|TH|) the factorization scale dependence of Tjj, H q , p o and <p p +. 

Eq. (|13|) describes the amplitude in the leading twist approximation. In other words 
all terms suppressed by powers of a hard scale parameter l/|p] are omitted. Within 
this approximation one neglects (in the physical gauge) the contributions of higher Fock 
states in the meson wave functions and many parton correlations (higher twist GPDs) 
in the proton. Moreover, we also use the collinear approximation: in the hard scattering 
amplitude we neglect the relative transverse momenta (with respect to the meson mo- 
mentum) of the constituent quarks. These approximations result in the appearance of 
the distribution amplitudes in the factorization formula (j!3|) . i.e. the lightcone wave func- 
tions, integrated over the relative transverse momenta of constituents up to the collinear 
factorization scale. 

An additional simplification appears in the kinematics given by Eqs. (J7HHJ)- In this 
limit one needs consider only the diagrams which involve two-gluon exchange between the 
two mesons. The other contributions (the fermion exchange diagrams) to the coefficient 
function are known [5] to be suppressed by powers of p 2 /s. Therefore we will not 
discuss them here. At the same accuracy, i.e. neglecting terms ~ p 2 / s, the contribution 
of gluon exchange diagrams shown in Fig. 2 turns out to be purely imaginary. It involves 
GPDs in the ERBL region — £ < x < £ only, which is quite specific for our process [4]. 

(*) 

The amplitude may then be written in terms of the impact factor J L / T L as 




i 




j^t^Pl (m £ ^ [H u (Z(2u - 1), e, 0) - H\i(2u - 1), £, 0)] . (14) 



o 



The impact factor J j l~*Pl has the form 



i 




(15) 



o 



with the 4>\\{z) DA defined by the matrix element 



i 




(16) 



o 



For 7W transversely polarized, J 7t ~* p ° l reads: 





(17) 
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Here e is the polarization vector of the initial photon, 
P(h,k 2 =p- h) = 9 ^ 9 — \-—^r- + 



z 2 p 2 + m 2 + Q 2 zz z 2 p 2 + m 2 + Q 2 ^ 
1 1 



(/ci — zp) 2 + m 2 + Q 2 zz (ki — zp ) 2 + m 2 + Q 2 zz 

and 

Qp(fci, h=p-h)= Y — — — — - -- ^—-JL— (19) 

z z p z + Q z z z + m z z z p z + Q z z + m z 

k\ — zp k\ — zp 



+ 



{ki — zp) 2 + Q 2 z z + m 2 (ki — zp) 2 + Q 2 z z + m 2 

In all our estimates in this paper we put the quark mass m q = 0. 

In the case of the transversely polarized pj meson production, one instead gets 



i 

X 





1 

/ dM f - a (n) rf^iuft up) [fl? (£(2t* ~ 1), 6 0) ~ iff (£(2tt - 1), e, 0)] , (20) 
J u z u z 



which involves the chiral-odd transversity distribution whose investigation is the main 
motivation of our studies. J^l/t^Pl are ^he same impact factors as in (|T3|) and 6 is the 
angle between the transverse polarization vector of the target n and the polarization 
vector e T of the produced p^-meson. In our numerical studies we take 9 = ir/2, but this 
dependence should of course be confirmed experimentally. 

The chiral-odd light-cone distribution amplitude for the transversely polarized p- 
meson is defined by the matrix element [6] 

i 

(p T (p p , T) | q(x)a^q(-x) | 0) = if J - [ due'^' 1 ^ . (21) 



We use the asymptotic forms of both pj and p° L DAs, i.e., <fi±(u) = 4>\\{u) = 6uu. 

We use the values of the meson decay constants in the above expressions at the scale 
1 GeV: /+ = 198 ± 7 MeV for the longitudinal p + , f° p = 216 ± 5 MeV for the longitudinal 
p°, and fj = 160 ± 10 MeV for the transverse p + . 

The generalized transversity distribution in the nucleon target described by the po- 
larization vector n M is defined by the formula [7] 

f d ±_ e ™P+z- {N(p2fj n )\q{- Z -) io + * q^)\N(p 2 , n)> = ^+1*0*, n)ia + V(p 2 , n)B*{x, £, t), 

(22) 
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Figure 2: The differential cross section for the photoproduction of longitudinally polarized p' 
and p + as a function of the squared transverse momentum ■pL for £ = 0.1, 0.3, and 0.5. 



where terms vanishing at t = t min (i.e. = 0) have been neglected. Our convention is 

Finally, two remarks are in order. Note that as \M. 1 2 , both from Eq. (|14j) and Eq. (J2(J|) . 
is proportional to s 2 , the s-dependence cancels in the full cross section (112)) . We want 
to stress that this is valid only for the high energy limit discussed in Sec. 2. If one is 
willing to study these processes at lower energy, one should include all polarization states 
of exchanged gluons 4 and the additional contributions related to quark exchanges. 



4 Longitudinally polarized meson case 

Let us first estimate the rate in the longitudinally polarized meson case. As in other 
p meson production processes, the dominant contribution comes from the unpolarized 
GPD and we will neglect the E(x,£,,t) contributions which vanish at t = t m i n . 

Our process selects the isovector part of this GPD. We use a standard description of this 
GPD in terms of double distributions, due to Radyushkin [8]. In general, it should be 
supplemented by a D-term [9] contribution. Since the chiral quark model estimates show 
that this latter term is almost flavor independent [10] we will not include it in our model 

4 At high energies the longitudinal polarization states of i-channel gluons (or in the Regge theory 
language the "nonsense" polarizations) give the dominant contribution and we omit the contribution of 
other polarization states. 
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Figure 3: The differential cross section for the photoproduction of longitudinally polarized p° 
and p + as a function of £ for p2, = 2, 4, and 6 GeV 2 . 



calculation involving isovector GPDs. We assume thus, 



H(x,tt) 



i + e j o 



dyF q - ■ -^ ,y,t 



(23) 



where the double distribution i 7 " 9 (X, Y, i) is given by the ansatz [8] 



F q (X,Y,t) 



m a(X)6 Y( 



l-X-Y) 



(24) 



In the expression (j24jl . q(X) is the quark parton distribution function, for which we use 
the parametrization of Ref. [11]. The t-dependence of the GPDs is given by the functions 
F q (t), which are related to the electromagnetic form factors of the proton and neutron, 
but since we are only evaluating the amplitudes at t = t min their form is not important. 

To compute the cross sections we must in general evaluate the integrals over u and 
z numerically; in the photoproduction case care must be taken to properly treat the 
apparent divergences at u = z and u + z = 1. 

The differential cross section for longitudinally polarized p + production is shown in 
Figs. El H3 and 0] for the real and virtual photon cases. We conclude from them that the 
photoproduction rate is much larger than the electroproduction rate (note that we did 
not include the additional suppression factor coming from the virtual photon flux and 
instead show the transverse and longitudinal photon cases separately). 
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Figure 4: The differential cross section for <y^^ T (Q)p — >• p\p^n for transverse virtual photon 
(left) and longitudinal virtual photon (right), plotted as a function of p T for £ = 0.3 and Q 2 = 
0, 1, 5, and 10 GeV 2 . 



5 Modeling the transversity GPD 

To estimate the rate of our process in the case of transversely polarized p + we need 
to formulate a model for the transversity dependent GPD HT{x,!;,t). From Eq. (}2T)j) 
it follows that we only need a model for the transversity GPD in the ERBL region, 
— £ < x < £. Unfortunately, even in the case of forward h\ structure function only very 
rough bag model estimates are available [12]. The bag model estimate of the leading 
transversity GPD was given recently in Ref. [13]. Recent progress was done also with 
lattice methods which have estimated its first x moments [14]. 

We propose a simple model inspired by Ref. [15] in which the E GPD was evaluated 
within the chiral approach, where a meson exchange dominates the physical process 
encoded in the GPD. The analogous meson pole approach was applied in the case of the 
forward transversity distribution in Ref. [16]. Below we generalize the pole model of the 
nucleon tensor charge developed in [16] to the non-forward kinematics. 

We start with the effective interaction Lagrangian 
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Na^d u A»N 



(25) 



in which gANN is the coupling constant determining the strength of the interaction of the 
axial meson A with the nucleon N. Inserting the interaction term i J d 4 x CanN into the 
S-matrix element in the left hand side of ()22|) and using the reduction formula one can 
separate the contribution to of the axial meson pole. In this way, for transversely po- 
larized nucleons described by the polarization vector St = (0, n±, 0), using the definition 
(I22|) of H T , we get the non-forward version of Eq. (7) of [16] as 



H T (x,0 



9ANNfA_ (A • Sj 

2M N m\ 



- 25 ) 



It 



(26) 



where A is the transverse part of the momentum transfer vector r (see Fig. and f^- 
is related to the A meson decay constant. We have here used the definition [6] of the 
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Figure 5: The transversity GPD Ht(x,£, 0) for various values of £ for x in the ERBL region, 
in our model (thicker curves) and in the model of Ref. [13] (thinner curves). 



transversely polarized axial vector meson distribution amplitude, 
(0\q(-l/2z)a a ^q(l/2z)\A(k,X)) = 



%f\ T [e a {\)kP -e p (\)k a ] / due i/2{l ~ 2u)k - z <Pi(u) 



(27) 



in which e(A) is the polarization vector of meson A with the momentum k, and 4>j_(z) is 



the distribution amplitude. In arriving at (j2T?j) it is useful to note that the definition 
can be rewritten using the Dirac structure a +J 75 instead of ia +l [7]. 

According to the model of Ref. [16] we now identify the scalar product (A ■ St) with 
the average of the intrinsic transverse momentum of the quarks: (A • St) — > l/2(k 2 L ). 
Also, the valence quantum number of t-channel isovector exchange leads to the iden- 
tification of the axial meson as A = (1235) . In this way, using the SU(2) relation 

(n\dOu\p) = (p\uOu\p) - (p\dOd\p), (28) 

we obtain as the final expression for the valence part H T = — H T (compare with 
Eq. (12) of [16]) 



H T (x,£,t) 



2^2M N ml 2£ 



(29) 
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Figure 6: The differential cross section for photoproduction of p° L , pj, as a function of "p\ for 
i = 0.1, 0.3, and 0.5. 



where, according to Eqs. (8, 9, 10) of [16], 

fl = —fa, , f ai = (0.19 ± 0.03)GeV 2 , 
5 

9b!NN = ^J^® aiNN ' 9aiNN = 7.49 ± 1.0 , (30) 

and 

(kl) = (0.58 + 1.0)GeV 2 . (31) 

(We use the lower value (k±) = 0.58 in our numerical estimates). In the following we 
assume that the distribution amplitude (j) b ^_ takes its asymptotic form, (j) b ^_ (u) = Quu. We 
show in Fig. El the resulting transversity GPD Ht{x,C,,0) for various values of £. 

The second model which we use in our estimates is the bag model estimate of Ref. [13]. 
In the Fig. we show also the comparison of this transversity GPD with the one defined 
by Eq. We see that these models lead to quite different results which consequently 
can serve as an estimate of theoretical uncertainties of our rate estimates. We want to 
stress that quark models, such as in Ref. [13], naturally underestimate the GPDs in the 
ERBL region, since they do not include the physics of meson exchange which is at the 
core of our model. 



6 Cross sections for p L pj, production 

Let us now estimate the rates for the production of transversely polarized meson and 
discuss the feasibility of measuring the transversity GPD in the proposed process. We 
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Figure 7: The differential cross section for the photoproduction of transversely polarized p { 
and p + as a function of £ for pL = 2, 4, and 6 GeV 2 . 



first note that since the analytic expression (J29J) for the chiral-odd GPD is rather simple 
we can compute, in the case of photoproduction, the integral over u in (l20|) analytically. 
One obtains 5 



\A~iv^PlP%n\ ■ ^16 7T 3 s a 2 e C F gwrg^g^j) I - j 1 

^ lQ-o = sm^ ^ ^-j ( 32 ) 



and therefore for the cross section 

da _ 729 ^a 4 s a em C 2 [g blNN fp ' fp&l)] 2 sin 2 # 



(33) 



dpldtdi N* Mlm\ x e(l + 0|pT ' 

For Q 2 > we evaluate the integrals numerically. The cross section for the transversely 
polarized p case has a characteristic sin 2 6 dependence (see Eq. (J20j) ). which we do not 
explore here. 

The differential cross sections for real and virtual photoproduction of transversely 
polarized p + , based on Eq. are shown in Figs. |H1 Eland [HI As for the longitudinal 
p\ case, we conclude from them that the photoproduction rate is much larger than the 
electroproduction rate (note again that we did not include the additional suppression 
factor coming from the virtual photon flux). Fig. H shows a comparison of the corre- 
sponding photoproduction cross sections computed using our model (J29|l f° r Ht and the 
bag model estimate of [13]. Fig. El illustrates the strong sensitivity of the production rate 

5 We emphasize that this result holds for our simple model estimate of the chiral-odd GPD, and not 
for a more complete treatment. 
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Figure 8: The differential cross section for ^, T [Q)p — > p^p^n for transverse virtual photon 
(left) and longitudinal virtual photon (right), plotted as a function of for £ = 0.3 and Q 2 = 
0, 1, 5, and 10 GeV 2 . 



with respect to the used transversity GPD, and also the different £ dependence of the 
two. 



7 Discussion and conclusions 

Up to now we have shown that the process discussed here has a sizable cross section 
with respect to the quite large luminosities of current and projected high energy electron 
accelerators, In addition, one should ask the question of the detection efficiency of exist- 
ing or planned experiments in the kinematical region where the two vector mesons are 
produced. A crucial point in favor of a good detection efficiency of our process is the fact 
that the two vector mesons are produced with large transverse momentum. Let us detail 
this statement. Detecting a quasi-forward diffractively produced p° is now a routine work 
for experimentalists at HERA and COMPASS, and we do not expect any problem on 
this point in a forthcoming electron-proton collider. The fact that we require that this 
p° has a finite transverse momentum of order 1 GeV or more will improve the detection 
efficiency Consequently, regarding the first (most rapid) vector meson, the two charged 
pions emerging from it should be well measured in modern detectors. 

The new aspect of the process that we have studied is the necessity to measure the 
charged p meson which is not diffractively produced and which travels in the vicinity of 
the proton. Since one needs to measure its polarization, a reasonable angular coverage of 
the two outgoing pions is required. Without performing a detailed study of the acceptance 
corrected cross section for a given experimental setup, we can however make some definite 
remarks. The fact that the p° meson has a sizable transverse momentum with respect 
to the quasi-real photon (i.e. lepton) beam, and the fact that the transversity GPD 
is peaked at small values of t favors the case of a p + meson with a similarly sizable 
transverse momentum (say greater than 1 GeV), which is shared between both emerging 
pions. This is good news with respect to the detection of these pions. To go further, 
one should distinguish between fixed target experiment such as Compass and collider 
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Figure 9: Comparison of photoproduction of pj, computed using our model for and the 
model of [13]. 

experiments such as HI or ZEUS. In the first case, the longitudinal momentum of the p + 
meson is (in the laboratory frame) of the order of ^Js^ ) which is a few GeV, so that the 
emission angles of the p + and hence those of the two n mesons are within the acceptance 
of the spectrometer. In the collider kinematics, the efficiency may not be as good since 
the vector meson is boosted along the direction of the proton beam. Note however that 
its longitudinal momentum is proportional to £, so that there is a small £ region where 
the emerging mesons should not be hidden in the dead cone around the proton beam. 
These positive remarks do not mean that a detailed analysis is not needed, neither that 
special care is not required to improve as much as possible the detection efficiency, but 
they allow us to be quite optimistic about the feasibility of the proposed experiment. 

We believe that the possible eRHIC machine may become the best place where the 
process discussed here could be measured, provided experimental setups allow a large 
angular coverage, ensuring a sufficient detection efficiency and a good control of exclusiv- 
ity. The JLab CLAS-12 upgrade probably will have good enough detection efficiency for 
observing the two rho mesons, but only for relatively low p T of the order of 1-1.5 GeV. 
Moreover the smaller energy available prevents the theoretical framework used here from 
being adequate and one needs to supplement our studies by adding contributions com- 
ing from other polarization states of exchanged gluons and the ones coming from quark 
exchanges. This we leave for a future work. 

The experimental measurement of the transversity GPD can thus supplement the 
intense present activity to unravel the transverse spin structure of the nucleon. An 
experimental determination of the transversity GPD Ht seems feasible in photo- or elec- 
troproduction at high energies if the accelerator luminosity is of the order of what is 
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anticipated at a future high-energy electron-proton collider. 

Acknowledgments 

We thank S. Scopetta for providing us with the parametrization of the GPD from [13]. We 
acknowledge useful discussions on experimental prospects with H. Avagyan, E. Burtin, J. 
Ciborowski, F. Kunne, A. Sandacz, L. Schoeffel and J. Ukleja. We also thank P. Hagler, 
D. Ivanov and O.V. Teryaev for helpful discussion. This work was supported in part by 
the Polish Grant 1 P03B 028 28, the French-Polish scientific agreement Polonium, the 
Joint Research Activity "Generalised Parton Distributions" of the European 13 program 
Hadronic Physics, contract RII3-CT-2004-506078, and the Swedish Research Council. 
L.Sz. is a Visiting Fellow of the Fonds National pour la Recherche Scientifique (Belgium). 

References 

[1] J. P. Ralston and D. E. Soper, Nucl. Phys. B 152, 109 (1979); 
X. Artru and M. Mekhfi, Z. Phys. C 45, 669 (1990); 
J. L. Cortes, B. Pire and J. P. Ralston, Z. Phys. C 55, 409 (1992); 
R. L. Jaffe and X. D. Ji, Phys. Rev. Lett. 67, 552 (1991). 

[2] V. Barone, A. Drago and P. G. Ratcliffe, Phys. Rept. 359, 1 (2002); 
M. Anselmino , |arXi v: hep-p h /0512140 

[3] A. V. Belitsky and D. Muller, Phys. Lett. B 417, 129 (1998); 
P. Hoodbhoy and X. D. Ji, Phys. Rev. D 58, 054006 (1998); 
M. Diehl, Eur. Phys. J. C 19, 485 (2001); 

M. Diehl, T. Gousset and B. Pire, Phys. Rev. D 59, 034023 (1999); 
J. C. Collins and M. Diehl, Phys. Rev. D 61, 114015 (2000); 
M. Diehl and P. Hagler, Eur. Phys. J. C 44, 87 (2005). 

[4] D. Y. Ivanov, B. Pire, L. Szymanowski and O. V. Teryaev, Phys. Lett. B 550, 65 
(2002) 

[5] L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); 

E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 44, 443 (1976); 
E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977); 
I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 

[6] P. Ball and V. M. Braun, Phys. Rev. D 54, 2182 (1996). 

[7] M. Diehl, Eur. Phys. J. C 19, 485 (2001). 

[8] A. V. Radyushkin, Phys. Rev. D 59, 014030 (1999). 

[9] M. V. Polyakov and C. Weiss, Phys. Rev. D 60, 114017 (1999). 



14 



[10] K. Goeke, M. V. Polyakov and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 
(2001). 

[11] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Phys. Lett. B 604, 
61 (2004). 

[12] R. L. Jaffe and X. D. Ji, Phys. Rev. D 43, 724 (1991). 

[13] S. Scopetta, Phys. Rev. D 72, 117502 (2005). See also B. Pasquini, M. Pincetti and 
S. Boffi, Phys. Rev. D 72, 094029 (2005). 

[14] M. GSckeler et al. [QCDSF Collaboration], Phys. Lett. B 627, 113 (2005). 

[15] L. Mankiewicz, G. Piller and A. Radyushkin, Eur. Phys. J. C 10, 307 (1999). 

[16] L. P. Gamberg and G. R. Goldstein, Phys. Rev. Lett. 87, 242001 (2001). 



15 



